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Employing first principles electronic structure calculations in conjunction with the frozen-magnon 
method we calculate exchange interactions, spin-wave dispersion, and spin-wave stiffness constants in 
inverse-Heusler-based spin gapless semiconductor (SGS) compounds Mn 2 CoAl, TLMnAl, Cr 2 ZnSi, 

TbCoSi and TLVAs. We find that their magnetic behavior is similar to the half-metallic ferro¬ 
magnetic full-Heusler alloys, i.e., the intersublattice exchange interactions play an essential role in 
the formation of the magnetic ground state and in determining the Curie temperature, T c . All 
compounds, except TLCoSi possess a ferrimagnetic ground state. Due to the finite energy gap in 
one spin channel, the exchange interactions decay sharply with the distance, and hence magnetism 
of these SGSs can be described considering only nearest and next-nearest neighbor exchange inter¬ 
actions. The calculated spin-wave dispersion curves are typical for ferrimagnets and ferromagnets. 

The spin-wave stiffness constants turn out to be larger than those of the elementary 3d-ferromagnets. 

Calculated exchange parameters are used as input to determine the temperature dependence of the 
magnetization and T c of the SGSs. We find that the T c of all compounds is much above the room 
temperature. The calculated magnetization curve for M 112 C 0 AI as well as the Curie temperature are 
in very good agreement with available experimental data. The present study is expected to pave the 
way for a deeper understanding of the magnetic properties of the inverse-Heusler-based SGSs and 
enhance the interest in these materials for application in spintronic and magnetoelectronic devices. 

PACS numbers: 75.47.Np, 75.50.Cc, 75.30 


I. INTRODUCTION 

Heusler compounds and alloys are a huge family of 
intermetallic compounds^ and they owe their name to 
the German metallurgist Friedrich Heusler who in 1904 
studied the thermodynamic properties of Cu2MnAl4 The 
development of the research fields of magnetoelectron¬ 
ics and spintronics^ intensified the interest in the half- 
metallic Heusler compounds—£ In half-metallic magnets 
a metallic majority-spin electronic band structure and a 
semiconducting minority-spin electronic band structure 
coexist^ Such compounds could lead to the creation of 
a fully spin-polarized current, maximizing the efficiency 
of spintronic devices— In addition, half-metallicity in 
Heusler compounds is always accompanied by the so- 
called Slater-Pauling rule where the total spin mag¬ 
netic moment scales linearly with the number of valence 
electrons in the unit cell— - — Although a large num¬ 
ber of half-metallic Heusler compounds compounds has 
been studied, still novel properties are discovered among 
Heusler alloys, e.g. ferromagnetic or ferrimagnetic semi¬ 
conducting behavior f i£ - — paving the way for diverse ap¬ 
plications in spintronics/magnetoelectronics^ 1 

A class of materials bridging the gap between half¬ 
metals and magnetic semiconductors are the so-called 
spin-gapless semiconductors (SGS); magnetic semicon¬ 
ductors where there is an almost zero-width energy gap 
at the Fermi level in the majority-spin direction and a 
usual energy gap in the other spin-direction. 22 - In SGS 
(i) the mobility of carriers is considerably larger than in 


usual semiconductors, (ii) excited carriers include both 
electrons and holes which can be 100% spin-polarized 
simultaneously, and (iii) a vanishing amount of energy 
is enough to excite majority-spin electrons from the va¬ 
lence to the conduction band. Several compounds have 
been identified as SGSj22r— and among them exist a few 
Heusler compounds ; 16 ! 32- — and it is M^CoAl, an in¬ 
verse full-Heusler compound, which has attracted most 
of the attention due to its successful growth in the form 
of films. First in 2008 Liu and collaborators synthesized 
using an arc-melting technique M^CoAl and found that 
it adopted the lattice structure of inverse full-Heuslers 
(see Fig. [T| for a schematic representation of the struc¬ 
ture) with a lattice constant of 5.8388 A and a total 
spin magnetic moment of 1.95 /zb per formula unit- 36 
The lattice of inverse Heuslers is known as the XA(or 
X a )-structure and it is similar to the well-known L2i 
structure of full Heusler compounds like Co 2 MnAl where 
only the sequence of the atoms in the unit cell changes. 
Moreover electronic structure calculations yielded a fer¬ 
rimagnetic state with a total spin magnetic moment of 
1.95 pi,B per formula unit and an antiparallel coupling 
between the Mn nearest-neighboring atoms— Although 
the calculated structure in Ref. HI is that of a SGS, au¬ 
thors do not mention it in their article— In 2011 Meinert 
and collaborators studied again theoretically this com¬ 
pound and almost reproduced the calculated results of 
Liu et al. using a different electronic structure method— 
They have also calculated the exchange constants show¬ 
ing that they are short range and the magnetic state is 
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stabilized mainly due to the direct interaction between 
nearest-neighbors and predicted a Curie temperature of 
890 K— But it was not until 2013, when Ouardi et 
al identified the SGS behavior of M^CoAl and have 
confirmed it experimentally in bulk-like pollycrystalline 
films.— The experimental lattice constant was found to 
be 5.798 A, the Curie temperature was measured to be 
about 720 K and the total spin magnetic moment per for¬ 
mula unit was found 2 /zb at a temperature of 5 K— Fol¬ 
lowing this research work, Jamer and collaborators have 
grown thin films of 70nm thickness on top of GaAs^ 
but these films were found to deviate from the SGS 
behavior— On the contrary, films -grown on top of a 
thermally oxidized Si substrate- were found to be SGS 
with a Curie temperature of 550 Kj - 41 First-principles 
calculations of Skaftouros et al identified among the in¬ 
verse Heusler compounds four additional potential SGS 
materials: T^CoSi, T^MnAl, T^VAs and C^ZnSi, the 
latter three being also fully-compensated ferrimagnets, 
and V 3 AI for which one V sublattice is not magnetic 
and the other two form a conventional antiferromagnet— 
The SGS character of T^MnAl was also confirmed by 
Jia et al — Wollman et al — confirmed the conclusion 
of Meinert et al. that direct exchange interactions are 
responsible for the magnetic order in M^CoAl study¬ 
ing a wide range of Mn 2 -based Heusler compounds and 
predicted a Curie temperature of 740 K using the spher¬ 
ical approximation , 44 Skaftouros et al. have discussed in 
detail the behavior of the total magnetic moment in in¬ 
verse Heusler compounds including the SGS materials * 14 
Galanakis and collaborators have shown that defects keep 
the half-metallic character of M^CoAl but destroy the 
SGS character— Finally, recent studies on the effect of 
doping of M^CoAl with Co, Cu, V and Ti^ as well as 
the anomalous Hall effect have appeared in literature— 

Motivated by these recent advances in SGS systems, 
the present work aims at a prediction and understanding 
of their magnetic properties at elevated temperatures. 
We employ density-functional theory at the ground state 
augmented by a Heisenberg model Hamiltonian for the 
prediction of the temperature dependent magnetization. 

We present calculations of exchange interactions, spin 
waves and temperature dependence of the magnetiza¬ 
tion in five inverse Heusler compounds known to present 
spin-gapless semiconducting behavior studied in Ref. (H: 
Mn 2 CoAl, Ti 2 CoSi, T^MnAl, T^VAs and C^ZnSi. We 
find that magnetic behavior of the SGSs is similar to 
the half-metallic ferromagnetic full-Heusler alloys, i.e., 
the intersublattice exchange interactions play an essen¬ 
tial role in the formation of the magnetic ground state 
and in determining the critical temperature, T c . Note, 
that even in the case of zero total spin magnetic mo¬ 
ment in the unit cell, the compounds under study are 
fully-compensated ferrimagnets and not conventional an- 
tiferromagnets, and thus the critical temperature should 
be called “Curie temperature” and not “Neel tempera¬ 
ture” . It turns out that the T c of all compounds is much 
above the room temperature. The calculated magnetiza- 



FIG. 1: Schematic representation of the lattice structure of in¬ 
verse Heusler compounds having the chemical formula X 2 YZ 
where X and Y are transition metal atoms (with the valence 
of Y larger than of X) and Z is an sp-element. On the right we 
present the nearest and next-nearest neighbors of an A site. 
Note that the large cube contains exactly four primitive unit 
cells. 

tion curve for M^CoAl as well as the Curie temperature 
are in very good agreement with available experimental 
data. The rest of the paper is organized as follows: In 
Sec. nn we present the computational method. In Sec. HU 
we present the computational results, and Sec. HU gives 
the conclusions. 


II. COMPUTATIONAL METHOD 
A. Crystal structure and ground state calculations 

Prior to discussing the structure of the compounds un¬ 
der study, we should note that M^CoAl has been already 
synthesized experimentally^ and for the other four com¬ 
pounds one should calculate the formation enthalpy to 
establish their possible experimental existence. In Refs. 
[I 3 I and [III authors studied the formation enthalpies for 
810 Heusler compounds and concluded that the forma¬ 
tion of the compounds is not favored with respect to the 
constituent elements only in the case where the heavier 
main-group metals such as thallium, lead and bismuth 
are involved. Thus we expect the compounds under study 
to be thermodynamically stable. 

All five compounds under study are called inverse 
Heuslers and crystallize in the so-called XA or X a struc¬ 
ture, the prototype of which is CuHg 2 Ti, but usually 
the sequence of the atoms follows the chemical formula 
X 2 YZ. X and Y are transition metal atoms with the va¬ 
lence of X being smaller than Y, and Z is a sp atom. 
A schematic representation of the structure is given in 
Fig. Q] There are four atoms along the diagonal of 
the cube following the sequence X-X-Y-Z and thus the 
two X atoms sit at sites of different symmetry. We 
use the superscripts A and B to distinguish them. As 
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shown also in the figure, each atom at the A or C site 
resides at the center of a cube with nearest neighbors 
four B sites and four D sites sitting at the corners of the 
cube; equivalently each atom at a B or D site has four 
atoms at A sites and four atoms at C sites as nearest 
neighbors. Concerning the next-nearest neighbors each 
A site has as second neighbors six C sites (and each C 
site has six A sites as second neighbors), and equiva¬ 
lent is the situation for the B and D sites. The envi¬ 
ronment of each site is important for the discussion of 
the magnetic properties of these compounds. We have 
used the lattice constants determined in Ref-HU using 
the full-potential nonorthogonal local-orbital minimum- 
basis band structure scheme (FPLO)i2 with the GGA 
exchange-correlation potential.— For Mn 2 CoAl the cal¬ 
culated value of 5.73 A is slightly smaller than the ex¬ 
perimentally determined lattice constants of 5.8388 and 
5.798 A in Refs.[H and[H, respectively. We compute the 
ground state electronic properties using the full-potential 
linearized augmented plane wave (FLAPW) method as 
implemented in the FLEUR code^i combined with the gen¬ 
eralized gradient approximation (GGA) to the exchange- 
correlation potential as parameterized by Perdew et al£2- 


Here M a is the integrated magnetic moment of sub¬ 
lattice a. In order to obtain the spin-stiffness D of 
the compounds we fit a linear or quadratic form Z?|q| 
or D|q| 2 respectively to the adiabatic magnon energies 
along a high symmetry line in the neighborhood of the 
T-point (in cubic systems, D is isotropic). A linear be¬ 
havior is present for the conventional antiferromagnets 
and we find it here for the compensated ferrimagnets. 
We should also note that within our formalism, where we 
consider an adiabatic approach for the magnons, we can¬ 
not study the Landau damping of the spin-waves induced 
by electron-hole excitations. In the compounds under 
study we expect the adiabatic approximation to provide 
reasonable results since there is a Stoner gap separating 
the magnon spectra from the continuum Stoner excita¬ 
tion spectra as in most half-metallic magnets (with the 
exception when the Fermi level is exactly at the higher 
energy edge of the minority-spin energy gap). 


C. Temperature dependence of the magnetization 
and T c 


B. Exchange constants and spin-wave dispersion 


Next we present the formalism on which the calcula¬ 
tions of the exchange constants are based. Our starting 
point is the classical Heisenberg model with unit vectors 
e nai pointing along the spin moments at positions R„ Q 
specified by a cell index n and magnetic lattice index a. 
The spin moments interact via exchange coupling param¬ 
eters and the exchange Hamiltonian H ex is 

H ex = — — 'y ' Jrnn e ma ' e n/3 ( 1 ) 

mnoi/3 


Note that the magnitude of the moments is incorporated 
in the parameters J. The exchange parameters are ex¬ 
tracted by a least square fit from ab-initio total energy 
calculations performed using the FLEUR code for a set of 
spin spirals with randomized wave- vectors— To reduce 
the computational cost we do non-self consistent calcula¬ 
tions and apply the approximation based on the magnetic 
force theore m 52 ! 53 to obtain the total energy differences 
from the differences in sums of eigenvalues. 

In order to obtain the adiabatic magnon 
dispersion^ - — we set up the spin-wave dynamical 
matrix A(q) at a Brillouin zone point q and solve for 
eigenvalues which provides the magnon frequencies: 


A a/3 (q) 



J«7(o)M 7 


J a7 (q)^ \ 

\Mp\\M a \ J 


( 2 ) 


J aP ( q) = E J 0 n COs[q ■ (Roa - R n/3 )] (3) 

n 


We employ the classical Monte Carlo technique to cal¬ 
culate the temperature dependence of the magnetization 
that is derived from the Heisenberg exchange Hamilto¬ 
nian ©. The technique provides an excellent estimation 
of the critical temperature, T) MC , from the position of the 
peak of the susceptibility as a function of temperature 
calculated as x(T) = [(M 2 (T)) — (M(T)) 2 ]//cbT, where 
(• • •) denotes thermal averaging over Monte Carlo steps, 
M(T) is the magnetization (before averaging), and ks is 
the Boltzmann constant. Use of a correction for finite- 
size effects, e.g. the fourth-order cumulant method^ 
can give a more accurate estimation of the critical tem¬ 
perature, but the correction here is small (of the order 
of a few K) since we are using large simulation super¬ 
cells (12 x 12 x 12 primitive cells corresponding to 5184 
magnetic atoms). We employ the Metropolis algorithm— 
and as a random number generator we use the Mersene 
Twister— The input comprises both the inter-sublattice 
and intra-sublattice exchange constants considering only 
the magnetic atoms, neglecting any contribution due to 
the interaction of the low-moment sp atoms. Moreover 
we ignore the contribution of the interstitial region— 

From the exchange constants the Curie critical tem¬ 
perature T c can be also estimated within the mean-field 
approximation (MFA). Actually, the critical temperature 
is given by the average value of the magnon energies 
which in MFA is the arithmetic average taking all the 
magnon values with equal weight. Thus is an arithmetic 
property that the MFA estimation is larger than exper¬ 
imental values^ - —. The MFA expression of the critical 
temperature for a multi-sublattice material like Heusler 
compounds has been provided in literature ] 61 ' 62 . 
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TABLE I: Calculated atom-resolved and total spin magnetic 
moments (in p b) for the five spin-gapless semiconducting in¬ 
verse Heusler compounds under study having the chemical 
formula X2YZ. The superscripts A and B distinguish the two 
inequivalent X atoms; we present the sum of the Z spin mo¬ 
ment and the interstitial spin magnetic moments; abs stands 
for the sum of the absolute atomic spin magnetic moments. 
Note that we have used the equilibrium lattice constants as 
calculated in Ref. [32]. 


Compound 

a(A) 

m IX A] 

m [ X B] 

rn m 

Z-\-inter ] 


-ftlfabs] 

M112C0AI 

5.73 

-1.52 

2.61 

0.98 

0.18 

2.0 

5.29 

TiaCoSi 

6.03 

1.41 

0.71 

0.39 

0.49 

3.0 

3.00 

TbMnAl 

6.24 

1.13 

1.00 

-2.59 

0.46 

0.0 

5.18 

Ti 2 VAs 

6.23 

1.04 

0.42 

-1.61 

0.15 

0.0 

3.22 

Cr 2 ZnSi 

5.85 

-1.59 

1.64 

0.03 

-0.08 

0.0 

3.34 


III. RESULTS AND DISCUSSION 

This section is divided into three parts. In the first part 
we discuss the ground state properties and magnetic mo¬ 
ments of the studied compounds. The second part deals 
with the exchange interactions, spin-wave dispersions, as 
well as spin stiffness constants. The last part focusses on 
the temperature dependence of the magnetization and 
critical temperature. Note that Ti 2 CoSi presents similar 
behavior to Mn 2 CoAl and Ti 2 VAs presents similar prop¬ 
erties to Ti 2 MnAl. Thus we will focus our discussion 
mainly on Mn 2 CoAl, Ti 2 MnAl, and Cr 2 ZnSi compounds. 

A. Spin gapless semiconducting behavior and 
magnetic moments 

The first step in the study of these materials is to es¬ 
tablish their ground state properties at 0 K as obtained 
from our first-principles calculations. In Fig. [2] we have 
plotted the density of states (DOS) projected on the tran¬ 
sition metal atoms for Mn 2 CoAl, Ti 2 MnAl and Cr 2 ZnSi 
and in Table U we present the atomic and total spin mag¬ 
netic moments for all five compounds. We do not present 
band structure plots since they are known from previ¬ 
ous works. 32 ’ 36 In all studied compounds we get a finite 
gap in the minority-spin band structure (negative val¬ 
ues of DOS) and a zero-width gap in the majority-spin 
band structure. In the case of Mn 2 CoAl and Cr 2 ZnSi 
the Fermi level is located at the middle of the minority- 
spin energy gap and for Ti 2 MnAl the Fermi level is lo¬ 
cated at the left edge of the gap. The position of the 
Fermi level within the gap is important with respect to 
the coupling between collective and single electron excita¬ 
tions discussed later. For all five compounds under study 
the total DOS per formula unit (f.u.) is similar to the 
ones calculated in Ref.fH) with a different full-potential 
method and thus spin-gapless semiconducting behavior 
of these compounds is a robust prediction of ab-initio 
electronic structure calculations. 

Among the five studied compounds, Mn 2 CoAl is a fer- 


rimagnet and Ti 2 CoSi a ferromagnet with total spin mag¬ 
netic moments per f.u. of 2 and 3 fj, b, respectively, and 
the other three compounds combine the SGS character to 
a fully-compensated ferrimagnetic state presenting zero 
total spin magnetic moments per f.u. as shown in Table 
HI In the case of the four ferrimagnetic compounds, the X 
atoms at the B sites couple : either (a) antiferromagneti- 
cally to the X atoms at the A sites and ferromagnetically 
to the Y atoms at the C sites (cases of Mn 2 CoAl and 
Cr 2 ZnSi), or (b) ferromagnetically to the X atoms at the 
A sites and antiferromagnetically to the Y atoms at the 
C sites (cases of Ti 2 MnAl and Ti 2 VAs). This behavior 
is expected from the so-called Bethe-Slater curved The 
early transition metal atoms like, Cr and Mn, when they 
are close to each other in space tend to have antiparal¬ 
lel spin magnetic moments. On the other hand for the 
nearest Ti atoms, the coupling tends to be ferromagnetic 
(Ti' 4 and Ti s atoms in Ti 2 YZ compounds). This behav¬ 
ior is also reflected on the exchange constants calculated 
and presented in the next subsection. Interestingly even 
in the case of Cr 2 ZnSi, where Zn is almost non-magnetic 
since all its d-states are occupied lying below the energy 
window presented in Fig. [2] the small induced spin mag¬ 
netic moment at the Zn and Si atoms leads to a small im¬ 
balance of the spin moments between the two Cr atoms 
(Cr' 4 has as a spin moment of —1.59 and Cr s of 1.64 
/is) and Cr A resolved DOS presented in Fig. [2] is not an 
exact mirror image of the Ci B DOS as in conventional 
antiferromagnets. 

Finally, we should shortly discuss the values of the 
atomic spin magnetic moments presented in Tablc[l] For 
all five compounds, spin moments are similar to the re¬ 
sults in Ref. [ 32 I where the FPLO(— electronic structure 
code has been employed. Moreover for Mn 2 CoAl results 
agree with the calculated values in Refs. I36B38I and for 
Ti 2 MnAl with the results presented in Ref. [42|. Since in 
each study a different full-potential ab-initio method has 
been used, we can be confident of the validity of our re¬ 
sults. Concerning now the experimentally available data, 
in Refs.[36| and [38| only the total spin magnetic moment 
per f.u. for Mn 2 CoAl has been measured which has been 
found to be 1.95 and exactly 2 /xb, respectively, in agree¬ 
ment with our calculated value of 2 ^b- A similar to¬ 
tal spin magnetic moment (1.94 fj, b) has been measured 
by Xu and collaborators^ for Mn 2 CoAl films on Si sub¬ 
strates at 5 K. The only discrepancy occurs in the study 
of Jamer et in Mn 2 CoAl films on GaAs(OOl), where 
X-ray magnetic circular dichroism (XMCD) experiments 
where carried out. Although XMCD is a powerful tech¬ 
nique, in Mn 2 CoAl there exist two Mn atoms with op¬ 
posite spin magnetic moments. XMCD can distinguish 
between different elements but cannot distinguish atoms 
of the same chemical element with different spin mag¬ 
netic moments and thus the values of the Mn moments in 
Ref. E3 cannot be interpreted as atomic moments. Note 
that in FLAPW method, the atomic spin magnetic mo¬ 
ments are calculated by integrating the spin-dependent 
charge density within each muffin-tin sphere surrounding 
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FIG. 2: Density of states (DOS) projected on the transition metal atoms. In the case of M112C0AI, positive (negative) 
DOS values correspond to the majority (minority) spin electrons. In the case of the other two compounds, which are fully- 
compensated ferrimagnets, spin-up (positive DOS) and spin-down (negative DOS) electrons have been chosen such that the 
atomic spin magnetic moments in Tablejl] have the right sign. Fermi level is set to the zero energy value. 


each atom. The interstitial region is not assigned to any 
atom. Although we have used almost touching muffin-tin 
spheres we can see in TablcQ]that hr the case of T^MnAl 
and T^CoSi a significant part of the spin magnetic mo¬ 
ment is located at the interstitial region. This should be 
attributed to the fact that early transition metals have re¬ 
states extending further away from the nucleus compared 
to late transition metals and the smaller the valence the 
more extended are the d-states. 


B. Exchange interactions and spin-wave dispersion 

Our calculations have shown that in the compounds 
under study inter-sublattice exchange interactions play a 
dominant role in formation of the magnetic state and 
critical temperature. In Fig. [3] we present the inter¬ 
sublattice exchange constants as a function of distance. 
Negative values of the exchange constants reflect an an¬ 
tiferromagnetic coupling of the corresponding spin mo¬ 
ments and positive values a ferromagnetic coupling. In 
all compounds the inter-sublattice nearest neighbor in¬ 
teractions dominate and especially the interaction be¬ 
tween the X s atom and its X A and Y nearest neighbors 
(Fig.Q]). The interactions between next-nearest neigh¬ 
bors X" 4 and Y are expectedly weaker. In the case of 
Cr 2 ZnSi the Zn atom is almost non magnetic while in 
the case of Mi^CoAl the intra-sublattice exchange con¬ 
stants between the Mn A -Mn 4 atoms have a sizeable value 
despite their large distance. In the case of T^MnAl, the 
Ti' 4,5 have positive spin moments and the Mn c has neg¬ 


ative spin moment as shown in Table[I] The situation is 
different to M^CoAl where Mir 4 site has negative, and 
Mn s site and Co c have positive moments. Thus the X" 4 - 
X s and X b -Y interactions have different signs for the two 
compounds. We can conclude that for the SGSs under 
study the interactions are short range: as seen from the 
Fig. [3] they decay quickly with the distance. This can 
be attributed to the existence of the finite spin-gap as 
discussed in literature in detail for half metallic Heusler 
compounds! 66 ' 67 

The short range nature of exchange interaction in 
M 112 C 0 AI and similar compounds has been shown by 
Meinert and collaborators in Ref.Hzj. Moreover, in the 
case of M 112 C 0 AI our results agree very well with the re¬ 
sults obtained in Ref.[37, considering the difference of a 
pre-factor 1/2 in the definitions of the Heisenberg Hamil¬ 
tonian between Ref.[!3 and the present work. (Also note 
that in their notations the sequence of the atoms is Co- 
Mn^-Mn^-Al and thus our Mn 4 atom corresponds to 
the Mn c atom in their article). Further calculations by 
Wollmann and collaborators on a series of Mn 2 -based in¬ 
verse Heusler compounds have confirmed the short-range 
nature of the interactions in these materials^ 

The contribution of each type of exchange interaction 
to the total exchange field Jo (the Weiss field acting on 
a spin moment) is given by the sum of the interactions 
over all possible pairs or, partially, over all pairs within a 
sublattice. The calculated values for the intra- and inter¬ 
sublattice Jo is presented in TablcHIl Our results confirm 
the conclusions drawn in the previous paragraph. The 
on-site inter-sublattice Jq are considerably larger than 
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FIG. 3: Inter-sublattice Heisenberg exchange parameters as a function of inter-atomic distances for (a) M^CoAl, (b) 
Ti2MnAl,and (c) Cr2ZnSi. 


TABLE II: On-site intra-sublattice ( e.g. J x ~ x = 52r. J^ R x , where R is the lattice vector) and inter-sublattice ( e.g. 
J x " x = R -x ) exchange constants (in meV), spin-wave stiffness constant D (in meV-A for compensated ferrimagnets 
or meV-A 2 for ferri- and ferromagnets), and the mean-field and Monte Carlo calculated critical temperatures (in K) for the 
compounds under study. Note that for the inverse X 2 YZ compounds the two X transition metal atoms occupy the A and B 
sites and the third transition metal atom Y occupies the C site. 


Compound 

jX A -X A 

Jo 

t x b -x b 

Jo 

7 Y-Y 

T x A -x B 

Jo 

tX A -Y 

Jo 

jX B -Y 

Jo 

D 

T c mfa (K) 

T c MC (K) 

Mn 2 CoAl 

-73.0 

-29.7 

-10.5 

-238.4 

-84.5 

132.4 

677 meV-A 2 

1134 

770 

TiaCoSi 

22.7 

6.9 

5.22 

169.4 

69.8 

0.6 

636 meV-A 2 

766 

550 

Ti2 M11AI 

6.5 

1.9 

-30.0 

101.4 

-184.4 

-196.8 

274 meV-A 

1222 

960 

Ti 2 VAs 

37.5 

2.3 

92.4 

45.8 

-145.5 

-55.1 

598 meV-A 

910 

800 

Cr 2 ZnSi 

21.2 

-18.1 

- 

-336.4 

- 

- 

752 meV-A 

1308 

1040 


the intra-sublattice Jo and only in the case of Mn 2 CoAl 
and Ti 2 VAs the jMn A -Mn A anc j jV-v ma kes a consider¬ 
able contribution into the total exchange field with dif¬ 
ferent signs. 

In Fig. [4] we present the spin-wave dispersion along the 
high symmetry lines in the Brillouin zone for Mn 2 CoAl, 
Ti 2 MnAl and Cr 2 ZnSi. Each spectrum has distinct fea¬ 
tures. First, the number of branches coincides with the 
number of magnetic atoms in the unit cell and thus we 
have three branches for Mn 2 CoAl and TI 2 MnAl and two 
branches for Cr 2 ZnSi where the Zn atom is not mag¬ 
netic. The energy dispersion curves of all compounds 
under study are typical for magnets with short range 
interactions, where nearest-neighbors and next-nearest- 
neighbors interactions dominate, and do not yield any in¬ 
stabilities. Instabilities can occur if the acoustic magnon 
modes have very low energies (close to zero) in some parts 
of the Brillouin zone but this is not the case for any of 


the studied compounds. 

In the case of Mn 2 CoAl the acoustic branch shows a 
typical behavior of ferro/ferrimagnets and around the 
r point the energy-dispersion curve shows a quadratic 
behavior with a spin-wave stiffness constant D of 677 
meV-A 2 . This value exceeds typical values of transition 
metal ferromagnets which usually range between 300 and 
600 meV-A 2 ,— and is close to the maximum known val¬ 
ues of 715 meV-A 2 for Co 2 FeSi^ and 800 meV-A 2 for 
Fe 53 Co 4 7 ^ Ti 2 CoSi as shown in Table HJ exhibits a D 
value of 636 meV-A 2 close to the value for Mn 2 CoAl. 

Although Cr 2 ZnSi is not a true antiferromagnet, 
around the T point the energy dispersion is linear and 
the optical and acoustic branch coincide showing consid¬ 
erable difference only at the plateau close to the Brillouin 
zone boundary (X, W, L, K points in Fig. [4]). The spin- 
wave stiffness constant as shown in TablclTD takes a value 
of 752 meV-A. In the case of Ti 2 MnAl the situation is 
more complicated. Due to the small spin moments of 
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FIG. 4: Spin-wave dispersion curves along the high symme¬ 
try lines in Brillouin zone for (a) Mn 2 CoAl, (b) Ti 2 MnAl,and 
(c) Cr 2 ZnSi. In the case of Mn 2 CoAl and TbMnAl, there 
are three branches since there are three magnetic atoms per 
unit cell; the sp-atom are almost non-magnetic. In the case 
of Cr 2 ZnSi, the Zn atoms have all their d- valence states com¬ 
pletely occupied and have vanishing spin magnetic moments, 
and thus there are two branches. In the case of Mn 2 CoAl the 
dispersion curve of the acoustic branch around the T point 
shows a quadratic behavior while for Cr 2 ZnSi the behavior is 
linear In the case of TbMnAl, it is in between. 


the Ti atoms around the T point the energy-dispersion 
of the acoustic magnon is practically linear but with a 
small value of the D constant (274 meV-A). T^VAs also 
shows a behavior similar to T^MnAl but now D has a 
much higher value of 598 meV-A. 


C. Temperature dependence of magnetization and 

T c 

In this section we show the temperature dependence of 
the magnetization and susceptibility calculated with the 
classical Monte-Carlo technique. In Fig. [5]we present for 
Mn 2 CoAl, Ti 2 MnAl and C^ZnSi the temperature de¬ 
pendence of the sublattice and total magnetization per 
f.u. First, we compare our results to experiment. We 
plot for Mn 2 CoAl our theoretical results together with 
the experimental results from the Ref.[38| where the tem¬ 
perature dependence of the total magnetization in a pol- 
lycrystalline film was measured in a field of 1 T. The 
agreement between the two data sets is good with one 
curve falling on top of the other with the exception of 
the region close to the critical temperature. Our curve 
shows an abrupt decrease close to our calculated critical 
temperature of 770 K, while in experiments this sharp de¬ 
crease is shifted lower in temperature since the measured 
critical temperature in Ref. [38| is 720 K, slightly smaller 
than our value. Thus the Monte-Carlo technique, which 
we employ, accurately describes the temperature depen¬ 
dence of the magnetization, indicating also the accurate 
determination of the exchange constants. 

In the case of C^ZnSi, the picture is that of an an- 
tiferromagnet with a total spin magnetic moment be¬ 
ing almost equal to zero for all temperatures, since we 
have ignored the magnetic properties of the Zn atom. In 
the case of T^MnAl the calculated total magnetization 
(Fig. [5]) is not exactly zero as expected from the first- 
principles result at 0 K (Table|l| but it equals to —0.46 
Hb due to the spin moment of the A1 atom and of the 
interstitial region which we have ignored in our Monte- 
Carlo calculations. At and close to room temperature 
(up to ss 400 K), where most devices operate, we find 
that all three compounds still have sizeable values of the 
sublattice magnetization. Thus the magnetic properties 
are not deteriorated although the spin-gapless behavior 
could be lost at this elevated temperature. For example 
for Mn 2 CoAl at 0 K the Mn A , Mn s and Co atoms show 
a sublattice magnetization of about —1.5, 2.6 and 1.0 /x b 
respectively. At room temperature, these values become 
— 1.3, 2.3 and 0.8 /xb, showing an absolute-value decrease 
of 13% , 11% and 18% , respectively. The total magne¬ 
tization per f.u. decreases from 2.0 /ib at 0 K to 1.8 /rs 
at room temperature showing an even smaller decrease 
of 10 %. Thus the compounds under study are adequate 
to be employed in spintronic/magnetoelectronic devices. 

To estimate the critical temperature T c MC within the 
Monte-Carlo technique we have plotted in Fig. [HI the sus¬ 
ceptibility x versus the temperature for all compounds 
under study. Obtained values are presented in Tablcllll 
and compared with the mean-field estimation of the crit¬ 
ical temperature T C MFA . As seen for all compounds T c MC 
is much lower than the T C MFA . The difference ranges 
from 110 K for T^VAs up to 364 K for M^CoAl, while 
for the other compounds it is around 220-270 K. Note 
that within MFA, the critical temperature is given by the 
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FIG. 5: Calculated temperature dependence of the sublattice and total magnetization for (a) M^CoAl, (b) T^MnAl, and (c) 
CrgZnSi. In the case of TbMnAl the total spin magnetic moment for 0 K is —0.46 p b since in the Monte-Carlo calculations 
we ignore the sp-atoms as well as the interstitial region (see Tabled]). With the arrow we denote the room temperature of 300 
K. In the case of Mn 2 CoAl we include also the experimental results from Re(.f38l. 



Temperature (K) 


FIG. 6: Calculated temperature dependence of the suscep¬ 
tibility for SGSs. The maximum of the susceptibility corre¬ 
sponds to the critical temperature (shown in parenthesis in 
the legends). 


arithmetic average of all the magnon energy values with 
equal weight. In reality low-energy magnons have more 
weight in determining the critical temperature and MFA 
usually overestimates experimental data by more than 
20 %j£& - — The Monte-Carlo determined critical temper¬ 
ature does not suffer from this drawback and thus it is 
expected to approach the experimental values of the crit¬ 
ical temperature. 

In some magnetic compounds, an empirical trend is 
observed that the critical temperature increases with 


the magnitude of the local magnetic momentsThis 
trend is not observed here. The sum of absolute local 
moments (see Table HI) is largest for M^CoAl being 5.29 
/is but the critical temperature is largest for C^ZnSi for 
which this sum is only 3.34 /tb ■ Intuitively we can explain 
this observation as follows. The critical temperature 
depends both (a) on the nature of the spin-dependent 
wavefunction overlap and chemical bonds between near¬ 
est neighbors, which is a short-range property and is in¬ 
dependent of the magnitude of the spin moments, as well 
as (b) the more distant interactions which are formed 
largely because of susceptibility effects (one atom po¬ 
larizes the Fermi surface electrons and this is felt by a 
distant atom through the electron propagation) which 
largely depend on the magnitude of the spin magnetic 
moments. Here, since the nearest-neighbor interactions 
are dominant and the long-range ones practically vanish, 
(a) applies but not (b). 

Finally, we compare our values for T c with experi¬ 
mental and theoretical results on M^CoAl compound. 
Our calculated MFA and Monte-Carlo values are 1134 
K and 770 K, respectively. Meinert and collaborators!! 
estimated theoretically a T C MFA = 890 K, smaller than 
our MFA result. Wollman and collaborators!! calculated 
jiMFA _ 935 py in-between our value and the value in 
Refill and within the more accurate spherical approx¬ 
imation (SPA)!! a value of 740 K, close to our Monte- 
Carlo value of 770 K.— Experimentally, for the bulk-like 
pollycrystalline film, Ouardi and collaborators measure 
for the T c a value of 720 Kj!! while the experiments by 
Xu and collaborators on thin films of M^CoAl gave a T c 













9 


value of about 550 Kj^i this discrepancy is expected since 
films present critical temperatures significantly smaller 
that the bulk samples. Thus our Monte-Carlo value of 
770 K is in good agreement both with the SPA calcu¬ 
lated value of 740 K of Wollmann et al 4^ and with the 
experimental value of 720 K measured by Ouardi et a7— 

IV. CONCLUSIONS 

We employed first principles electronic structure cal¬ 
culations in conjunction with the frozen-magnon method 
to calculate exchange interactions, spin-wave disper¬ 
sion, and spin-wave stiffness constants in inverse-Heusler- 
based spin gapless semiconductor (SGS) compounds 
Mn 2 CoAl, T^MnAl, C^ZnSi, T^CoSi and T^VAs. We 
find that the magnetic behavior of the SGSs is similar to 
the half-metallic ferromagnetic full-Heusler alloys, i.e., 
the inter-sublattice exchange interactions play an essen¬ 
tial role in the formation of the magnetic ground state 
and in determining the critical temperature, T c . All com¬ 
pounds, except Ti 2 CoSi, possess a ferrimagnetic ground 
state. Due to the finite energy gap in one spin channel, 
the exchange interactions decay rapidly, and hence mag¬ 
netism of these SGSs can be described considering only 
nearest and next-nearest neighbor exchange interactions. 
The calculated spin-wave dispersion curves are typical 
for ferrimagnets and ferromagnets, i.e., the number of 
inequivalent spin-wave branches in the dispersion curves 
is equal to the number of magnetic atoms in the unit cell. 
Due to the short range nature of the exchange interac¬ 


tions, the calculated spin-wave stiffness constants turn 
out to be larger than the elementary 3d-ferromagnets. 

Calculated exchange parameters are used as input to 
determine the temperature dependence of the magneti¬ 
zation and the critical temperature T c of the SGSs. We 
find that the T c of all compounds is much above the 
room temperature. The calculated magnetization curve 
for Mn 2 CoAl as well as the critical temperature are in 
good agreement with available experimental data. 

Our results suggest that, except M^CoAl which has 
been already synthesized, there are other potential SGS 
presenting also very high Curie temperatures. In SGS 
materials, only one spin channel contributes to the trans¬ 
port properties, whereas the other spin channel allows for 
tunable charge-carrier concentrations. Among the five 
studied SGS compounds Ti 2 MnAl, Ti 2 VSi and Cr 2 ZnSi 
present also zero magnetization thus creating vanish¬ 
ing stray fields leading to advantages that have been 
pointed out for antiferromagnetic spintronic devices! 74 ’ 75 
Thus we expect such fully-compensated ferrimagnetic 
compounds to be more adequate for room temperature 
spintronic/magnetoelectronic applications based on spin 
transport. 
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